Inter-hemispheric asymmetry of echoic memory traces with multiple lifetimes 
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1 Introduction 

Earlier studies [1,2] suggested exponentially 
decaying memory trace with significantly higher 
lifetimes for the association cortex than the primary 
cortex, and little or no inter-hemispheric 
differences. In these studies only the peak MEG 
signal (where signal to noise ratio is largest) or 
peak strength of the equivalent current dipole was 
used to model the decay of the echoic memory and 
quantify its dependence on the inter-stimulus 
interval (ISI). A unitary cortical response to a 
stimulus was implied since the shape of the 
neuromagnetic waveforms was assumed to be 
insensitive to the ISI. Under this model, the 
physiological memory trace once established 
decays passively and uniformly in each cortical 
area with time thereafter, justifying the use of 
separate but latency independent lifetimes for the 
echoic memory in the primary and association 
auditory cortices. 

A more general approach was recently proposed 
[3], which computes the fit to the exponentially 
saturating function over the wider latency domain. 
This allows assessing latency-dependent effects that 
may involve different neuronal pools with different 
onset latencies and/or memory trace lifetimes. The 
present study refines the earlier analysis revealing 
new insights about the latency and especially 
interhemispheric dependence. The auditory stimuli 
sequence was modified so that the average response 
following each ISI comes out from the same 
number of single trials and mixes together the 
different past ISIs histories with equal 
contributions. Our new protocol provides a more 
objective evaluation of the evoked responses and 
better dissociates the echoic memory trace decay 
from refractory effects, which may easily dominate 
runs with constant ISI. As in one of our previous 
work [4], the present study makes use of a virtual 
sensor computation, which maintains the analysis 
tractable while saving computational effort. 


2 Methods 

2.1 Data acquisition 

The CTF whole head system was used to record 
MEG signals from 7 subjects (6 right handed and 
one left handed) with no history of otological or 
neurological disorders and normal audiological 
status (air conduction threshold no more than 10 dB 
hearing level). The subjects were instructed to 
listen passively to binaurally presented auditory 
stimuli arranged in four runs. The stimuli were 
delivered using echo-free plastic tubes. Each run 
contained 30 repetitions of a block of 9 tones. In 
each block the 9 tones were spaced in a fixed 
sequence using 6 ISI values (1, 1.5, 2, 3, 5, 8 sec). 
Each tone (1 kHz, 60 dB nHL) had a duration of 
200 msec (rising and falling times of 10 msec). 
Figure 1 shows one such block. A heavy arrow is 
printed after each 5 sec ISI (relative to the directly 
preceding tone) and the last stimulus in the block 
(which is also the first stimulus of the next block) is 
distinguished with a dashed outline. The sequences 
in the blocks were designed to allow the cross¬ 
averaging of the responses following the same 
immediately (first) preceding ISI but with 
equiprobable presentation of a different history of 
the past (second preceding) ISI. The short duration 
of each run and the passive listening paradigm were 
used to avoid learning of the block-sequences by 

the subjects. We denote by s x \ y the evoked 

response elicited by a tone stimulus, which belongs 
to the nth block of a run, it is preceded by an inter¬ 
stimulus interval x and conditioned by a past 
(second preceeding) inter-stimulus interval y. This 
allows us to define the intra-run sub-averaged 
response s x \ y as: 
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*x\y = T7 M|v > N=30. (1) 
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The cross-averaged (inter-run) response to a first 
preceding ISI x is simply defined by: 
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Figure 1. The design of the ISI sequences. Each ISI 
value is preceded once and only once by each ISI 
value in the set, as shown by the s^ y example (note 

the double occurrence of s^ y in the first and the 
third sequence). 

The design of the ISI sequences ensures that 
different histories have equal contributions to the 
cross-averaged signals (each ISI is once and only 
once preceded by each of the ISIs values from the 
list, as shown in figure 1 for evoked responses 

)• 

2.2 Virtual sensor computation 

The early peaks in the average MEG signal elicited 
by auditory stimuli, particularly the Ml00, have 
predominantly dipolar distributions. A weighted 
average can be defined around each channel, i.e. for 
channel i: 

VS t (t)= £ exp[-(T-fk) 2 ]. Sj (t), (3) 

j=l A 

where X is the characterstic length and Sj is the 
MEG signal measured at channel j, whose position 
vector is r 7 . This smoothing operation minimises 
noisy contributions and can be applied separately to 
the extrema of the dipolar pattern [4]. Even more 
simply the five channels (ki to k 5 ) that produced the 
most clear positive deflections and five channels (pi 
to p 5 ) that produced the most clear negative 
deflections at the time of the Ml00 peak can be 
selected to define the composite virtual sensor 
signal (henceforth refered to as VS): 

VS(t) = -[ £ VS k . (t)-'tVS P . (t)J . (4) 

•> i=l i=l 



Figure 2. Example of channel selection for VS 
computation. Dashed lines surround the selected 
channels with positive deflections, while continuous 
lines surround the selected channels with negative 
deflections. Lower panels illustrate the resulting 
left and right VS signals. 


We have adopted this model-independent procedure 
to define one pair of VS for each subject in each 
hemisphere. The effectiveness of the operation is 
exemplified in figure 2. In our earlier studies [3,4] 
we have shown that the virtual signal is dominated 
by contributions from generators associated with 
activity in the auditory cortex. This was achieved 
by moving each 37 channel probe of the twin 
MAGNES BTI system over the expected auditory 
cortex position, separately in each hemisphere, until 
the extrema of the dipolar pattern were clearly 
captured. In the new experiments we used the 151 
whole head CTF system which provides even 
coverage over the head. The sensors are further 
removed from the brain and hence are less spatially 
selective than in the earlier experiments; the 
selection around the Ml00 and the overall location 
biases the VS in this study in favour of activity 
generated in the auditory cortex. Although possible 
to design more selective combinations and target 
the auditory cortex directly, we chose to confine 
ourselves to the model independent approach we 
had previously used. 

2.3 Data processing 

The continuous data stream in each VS was first 
divided into 700 msec segments, each beginning 
200msec before stimulus onset. The intra-run and 
inter-run averaging was then performed as 
described above. In this way, different past histories 






















































































Figure 3. The cross-averaged VS signals for subject 2 (NL), left auditory cortex (a), the corresponding fitting 
parameter z versus latency (b), and the curve fitting at latencies 112 msec (around M100) and 136 msec (c). 
The cross-averaged VS signals for subject NL, right auditory cortex (d), the corresponding fitting parameter z 
versus latency (e), and the curve fitting at latencies 112 msec and 136 msec (f). 


were mixed together and the dependence of the 
response to only the first preceding stimulus is 
emphasised. A 4 sample moving average was used 
to smooth out the small high frequency 
fluctuations. The signal was further processed at 
latencies with significant and consistent activation, 
selected on the basis of the following simple rule: 
the signals from all ISI must be either all positive or 
all negative (intuitively, a negative value in a set of 
positive values for example will provide an 
unacceptable fitting). At these latencies, the least 
square best-fit estimate has been computed for the 
three parameters of the model function: 

t 0 -ISI 

A(1 - e T ) (5) 

where A is the amplitude, t 0 is the intercept or time 
of decay onset, and z is the lifetime. The fitting 
criteria were further restricted such as the 
convergence of the fitting algorithm to a tolerance 
of 0.001 is achieved within at most 10 iterations. 
The goodness-of-fit has also been visually assessed 
at each latency that past the previous criteria. 

In our earlier study [3] we have allowed small 
lifetime values, for which our recent work has 


revealed unstable fits arising from flat distributions. 
We have therefore excluded lifetime values smaller 
than 0.5 sec, for which a separate analysis effort is 
ongoing using sets with ISIs less than Is. 

3 Results 

Figure 3, panels a and d illustrate the cross- 
averaged VS signals for subject 2, left and right 
auditory cortex, respectively. Around each major 
positive and negative deflection of the VS signal 
the dependence of the activation strength on ISI is 
evident, and it is well described by the 
exponentially saturating function. Panels b and e 
show that the latency dependence of the parameter 
T splits into segments that repeat around each peak 
of the signal. In the left hemisphere lifetime values 
in each segment increase with latency, while on the 
right a segment has an inverted “U” shape. Panels c 
and / shows the VS data and the corresponding 
fitting function at two latencies, i.e. around the 
Ml00 peak (latency 112 msec) and on the 
descending part of Ml00 (latency 136 msec). The 
lifetime values are almost the same for the two 
hemispheres around the peaks (1.61 and 1.36 sec). 








































Table 1: The T values at the mean Ml 00 latency 
(Lmioo) an d the maximum T values in the range 
100 to L m100 +30 msec. The first subject is 
the only left-handed (female) member of the set. 



T values (sec) at 
latency L M100 

Maximal T values 
(sec) between 

L M100 and 

L m joo +30 msec 

Subject 

Left 

Right 

Left 

Right 

i 

1.12 

1.31 

3.68 

1.85 

2 

1.61 

1.36 

5.45 

2.24 

3 

0.95 

0.7 

6.01 

1.11 

4 

0.92 

0.84 

5.24 

3.64 

5* 

0.88 

0.61 

4.47 

1.01 

6* 

1.46 

1.15 

4.05 

1.21 

T 

2.41 

1.82 

4.69 

2.91 

Mean±SD 

1.3±0.5 

1.1±0.4 

4.8±0.7 

1.9±0.9 


Only on the left hemisphere T assumes significantly 
higher values at latencies immediately following 
the peaks, especially the Ml00. The intercept t 0 has 
been found to have different values with latency, 
which can no more be associated with the stimulus 
duration as suggested in [1]. Table 1 shows the 
consistency of the inter-hemispheric differences in 
lifetime values across subjects, by contrasting the T 
values corresponding to the left and right auditory 
cortex at the mean latency of the Ml00 peak 
(denoted by L M jqq) and the maximal values in a 
time window between L m10 q and L m 1 qq+ 30 
msec. For some subjects (marked with a star in 
Table 1), the signal for ISI=8 sec was very low; 
while this phenomenon remains unexplained, the 
corresponding data have been considered as deviant 
and removed from the analysis before passing the 
data to the fitting routines. 

4 Discussion 

The present study confirms that the morphology of 
the auditory evoked response (in contrary to only 
the peaks) may encode information about the 
history of preceding activations. The encoding 
factors can be thought of as memory traces with 
different lifetimes. Our findings are consistent with 
the hypothesis that memory traces encoding 
physical features of the acoustic input are 
associated with the generation of the mismatch 
negativity (MMN), rather than neurons sensitive to 
the changed stimulus. 


The inter-hemispheric difference we have identified 
can arise from a single neuronal pool, sequentially 
modulated by a complex mechanism, which 
integrates memory traces with different lifetimes. 
Alternatively it may be generated by sequential 
adjacent sources characterized by memory traces 
with different characteristic lifetimes. The model of 
sequentially activated adjacent generators (partially 
temporally overlapping) of the Ml00 complex has 
already been proposed in [5] and its inter- 
hemispheric asymmetry has been suggested in [6]. 
Inter-hemispheric asymmetry is not evident at the 
signal peaks, and this explains why previous studies 
failed to recognised it. In our study it is identified 
for the first time and it is most evident at latencies 
immediately following a peak. The left hemisphere 
is characterized by significantly higher lifetimes of 
the memory traces so that it is more appropriate to 
integrate complex information in a sequential 
fashion. The left hemisphere was also higher for the 
only left-handed subject in the study, but with the 
lowest maximum lifetime across subjects on the 
left. It is not clear at this point what if any are the 
functional correlates of dependence of the lifetimes 
on latency, or if and how lifetime values change 
with age and in pathological. 
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